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Abstract 

We propose a new, very flexible version of the Rayleigh-Schrodinger perturbation 
method which admits a lower triangular matrix in place of the usual diagonal un- 
perturbed propagator. The technique and its enhanced efficiency are illustrated 
on rational anharmonicities V^^\x) = P x polynomial(x)/polynomial(x). They are 
shown tractable, in the intermediate coupling regime, as 0{P — P^^^) perturbations 
of exact states at non-vanishing P^^^ ^ 0. In this sense our method bridges the gap 
between the current weak- and strong-coupling expansions. 

PACS 03.65.Fd 03.65.Ge 02.60.+y 02.70.+d 



1 Introduction 

Schrodinger equations with anharmonic potentials uP'x^ + /3 V^^^ [x] are often solved 
perturbatively. It is well known that many practical implementations of this ap- 
proach are full of contradictions, well illustrated by the popular quartic example with 
y(i)(x) = x^. Its weak-coupling energy estimates E{(5) ^ E(0)+/5-E(i) + . . .+/9^-E(Ar) 
are easily generated via recursion relations [^ but this approximation diverges in the 
limit A^ ^ cxD at any nonzero coupling constant /5 . An alternative, strong-coupling 
series in powers of /5~^/^ exists and converges for the sufficiently large |/3| 0. Un- 
fortunately, the explicit evaluation of its coefficients is by far not easy 0]. In the 
literature many people have advocated, therefore, a replacement of the traditional 
quartic model by a non-polynomial anharmonicity 

Its merits belong to the two separate categories. Firstly, its bounded character 
enables us to avoid the divergence of the weak-coupling series. This has been em- 
phasized by several authors . Secondly, the existence of a few elementary solutions 
at certain exceptional couplings /3 = P^^^ [Q enables us to contemplate their pertur- 
bations, say, 

E{I3) = E{I3^''^) + AE(i) + A'E^^) + . . . , X = (3- /?W (2) 

in anharmonic regime, near any solvable P^^'' 7^ [0. An exceptional character of the 
latter intermediate-coupling expansion was its shortcoming. Its feasibility relied on a 
replacement of the traditional unperturbed spectrum by certain auxiliary continued 
fractions and did not seem amenable to any sufficiently efficient generalization [^. 
In the present paper we shall consider the whole class of the Pade-like forces 

V{(3, x) = u'x' + (3 (J2Bd xA E A^ x^", B, ^ 0. (3) 

\d=0 / n=0 

Their strongly anharmonic perturbation solutions of the type (j^) will be described 
in full detail. We shall demonstrate that near many non-vanishing "intermediate" 



couplings P^^^ our perturbative approximants may be both feasible and quickly con- 
vergent. 

The project opens several technical questions. Firstly, at any t > 1, the zero- 
order solutions must be constructed in anharmonic regime. A representative sample 
of these reference systems is described thoroughly in Section 2. It underlines the real 
phenomenological appeal of eq. (^) in comparison with the more popular polynomial 
models. 

At a particular anharmonic /3''°^ 7^ the usual construction of a complete unper- 
turbed basis is prohibitively complicated. After a return to harmonic basis, numer- 
ical integration is needed for evaluation of the necessary matrix elements of V^^^ (x) 
and the computation of corrections is difficult even in the lowest order Rayleigh- 
Schrodinger approximation 0. Section 3 offers the remedy. Schrodinger equation is 
represented in a non-orthogonal basis. Its resulting (2t -|- 1)— diagonal matrix form 
is then much better accessible to a purely numerical matrix-inversion perturbative 
treatment. 

In our main Section 4 we deny the latter numerical "brute force" philosophy and 
intend to soften it significantly. Re-installing the more traditional recurrent interpre- 
tation of perturbation algorithms we describe a new approach to the Schrodinger-like 
families of equations with a banded-matrix form of their Hamiltonians. Our main 
idea is amply illustrated by its application to anharmonicities (|1|). Its core is a maxi- 
mal simplification of the unperturbed propagator TZ. In contrast to its general- matrix 
form in older methods ||10| we shall be able to reduce it to the mere "half- filled" , 
triangular matrix. 

Section 5 is the summary showing how our new approach opens a way towards a 
broader variability of shapes of the theoretical and/or phenomenological interaction 
models. With due attention paid to the non-hermiticity of our (quasi-)Hamiltonian 
matrices (cf. also Appendix), our new version of perturbation recipe seems well 
prepared for its extensions as well as further practical computational applications. 



2 Solvable oscillators with /3^o) 7^ 

2.1 The simplest model with t = 1 

Non-polynomial eq. (|^) with t = 1 is often recalled as one of the simplest unsolvable 



anharmonic models in one dimension |rT|. In place of using the differential form 
of its Schrodinger equation, wave functions ^(/3,a;) G L2(— 00,00) are expanded in 
the harmonic (i.e., Hermite or Laguerre) polynomial basis {\n)}'^^Q. The ansatz 
ipi^P, x) = {1 + B x^) J2'^=o {x\n) hn and scaling cu ^ 1 then give ||12[ the three-term 
recurrences 



/ 



ao do 

C]^ di d]^ 

C2 02 



V 





( hA 




hi 




h2 




{ ■■ J 



(4) 



an = (3 + {1 + Ban) (Sn -E), an = en/2 = (n|r» = 2n + i + 3/2, 
dn = B(3n (sn ' E), /3„ = (n|r> + l) = [(n + l)in + i + 3/2)]i/2 

Cn = B Pn-l{£n- E), n = 0, 1, . . . . 

Marginally, let us note that the parity (—1)^"^^ of the wave functions with i = 



-1, 



admits an immediate re- interpretation as angular momentum i = 0, 1, . . . in three 
dimensions with the regularity 'ijj{(3,r) ~ r^"*"^ of the radial wave functions near the 
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2.1.1 Termination conditions and exact solutions 

The existence of the terminating exact solutions of eq. (^ is well known 0. Let 
us mark them by a superscript ^^'. With normalization h''^' ^ their termination 
property 

•• = (5) 



/.(O) _,(o) 

''g+l ~ "'13+2 



requires that c^Iy = 0. This means that E^'^^ = Aq + 2£ + 7 = e^+i. The related 
energy is not arbitrary. Vice versa, bound states ip{P,x) with E ^ E^^^ have to 
be defined by infinite series |T^. In one dimension a sample of their spectrum is 



given in Table 1 . It indicates that with a growth of the barrier the low-lying energies 
merge in almost degenerate doublets with opposite parities. According to Figure 1 
a very good fit of these numerical values E = E{(3) is provided by parabolas. One 
may expect that besides our ansatz (^ a useful methodical alternative could be also 
sought in perturbative expansions of couplings (3 = I3{E) and of the related Sturmian 
wave functions [l^, |15[. 

Even for the terminating bound states with E = E^^^ we have to guarantee that 
the secular determinants vanish. Up to g = 3 the latter condition is non-numerical. 
For illustration we may fix B = S^'^^ = 1 and choose the even parity i = —1. Then 
we get the elementary implicit polynomial definitions 



(6) 



y-6 = 0, q = 

?/2 - 26 y + 152 = 0, g=l 

y3 - 68 y2 + 1372 y- 8304, q = 2 

y^ - 140 y^ + 6588 y^ - 123216 y + 777600 = 0, g = 3 

of the partially solvable couplings y = y{q) = P^^\ Besides their g = (linear), g = 1 
(quadratic) and g = 2 (Cardano) explicit solutions we may write down all the four 
exact g = 3 roots 



/^W = 35 + £1 Jl27 - 2 V6821 cos (-Q + -] + £2 Jl27 + 2 V6821 sin (-Q + - 



+£3 1/127-2^6821 sin-G, 9 = arctg 



9 V63908723442661575155 



3 ° 45902084710 

with (ei, £2, £3) = (-, -, +), (-, +, -), (+, -, -) and (+, +, +). 

All our termination-compatible g < 3 values of the coupling /5*^°^ remain real (cf. 
their list in Table 2). All the related forces acquire a double well shape since all our 
roots satisfy its sufficient condition /?(°) > 1. 

2.2 Unperturbed solutions with t = 2 

The t = 2 option in (|^) gives the "first nontrivial" potential 
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V{x) =x^+ -\- ^ . (7) 



Mutatis mutandis, equation (H) becomes replaced by the g + 3 relations 
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The last row implies that d^'^2 — ^^d fixes the energy i?*^°^ = 4g + 2£ + 11. The 
coupled rest remains over-determinate and defines the g + 1 unknown coefficients 
(normalized, say, to /i^°^ = 1) and two coupling constants. At g = we have 



V 



(0) 



(4£ + 6) (/ - 2(?) + 8, /i^"^ = (8£ + 20) <7' + 4/ - S^? g = . (9) 



These couplings are real for all the parameters /, g and i. 

For the sake of brevity, let us put f = g = 1 from now on. With option i = 
(in one or three dimensions) and degree g = 1 our conditions (§) degenerate to the 
cubic equation u^ + ASu — 360 = 0. Its only real Cardano root 
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z/(°) = (4V228I + I8O)''''- (4^2281-180]"" ^4.95914661133166 (10) 

with /i^°) ^ 14.941997536546 and normalization h\ = 1 leads to the exact 
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i.e., hf) ^ —3.48195017221496. The related energy E^^^ = 15 corresponds to the first 
excited state in s— wave. Its ground state predecessor does not terminate, g —>■ 00. 
For it, the Runge-Kutta integration gives Egs ~ 10.943408413. 

At the next choice of g = 2 with the quasi-harmonic energy E'^'^^ = 19 and with 



the same convenient normalization /ig = 1 equation (P) is fully pentadiagonal, 

/^|/i + i/-52 ^/i-32V6 -8^30 ^ 

^/i-24V6 |/i + z/-195 v^/i-96V5 

^5/2-64^5 ^/i + z/-330 

-2 V2T0 ^ /i - 24 V42 y 

In an ascending order of its rows we eliminate 
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With hi = y v5 our problem degenerates to the single sextic polynomial equation 
with integer coefficients, 

2500 / + 1000 / - 7125 y^ + 100 y^ + 5065 / -26Ay- 1152 = 0. 

Its two real roots are easily localized numerically, yi ~ —1.69830 and ^2 ~ 0.815078. 
A close analogy with t = 1 is preserved. With the same ease we may generate the 
two t = 2 oscillators from a tenth-degree polynomial at g = 3 (see Table 3) etc. The 
related energies were evaluated numerically. Their sample is given in Table 4. They 
safely stabilize at cut-off M = 15. 



2.3 t = 3 and more 

At any t > 1 the requirement (^) leads to the (t + g + 1) x (g + 1)— dimensional 
generalization of eq. (^. Its last, decoupled condition c[\^f = is satisfied if and 
only a E^^^ = 4t + 4g + 2£ + 3. The remaining t + q coupled equations 
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determine all the q normalized projections M plus t parameters in potential itself. 
With the growth of t the selfconsistent search for these exact solutions becomes less 



and less straightforward. Due to the implicit nonlinearity of eq. (|T3D we must verify 
that its solutions keep the forces real and non-singular. Both these properties have 
to be verified a posteriori. 

It is useful to notice that at g = the explicit solutions remain elementary at any 
index t > 1. At t = 3 the purely non-numerical solutions still exist at g > 0. This is 
slightly unexpected. For illustration, let us employ the quartic-over-sextic model 

^(/5(o) , r) = a;2 + ^^— ^ . (14) 

In a search for its symmetric bound states in one dimension (£ = —1) the "first 
nontrivial" choice of g = 1 gives E^^^ = 17. The abbreviation Kq = a and the 
eliminations guided by our previous experience re-parametrize the couplings, 
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2(^2 + 27), t;(°) = -72(72 0^ + 110-6 72) 



u 
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72 (2 a^ + 10 72 a^ - 23 a + 18 72)/2. 



The whole algebra degenerates to the single equation in y = 72 a, 

y^ + 9y^-33y2 + 27i/-12 = 0, £=-1, g = l, t = 3. 

Two of its roots are complex, 1/34 = 0.451 ± 0.534 z, and the real doublet is given by 
the expression 
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The last item is a positive square root of another sum. 
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'^'> = ^r + V ws + V ws 
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+V m + V m 
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^V 1024 ^ V 1024 ■ 

In the second illustration with q = 2, E^^^ =21 and abbreviations h^ = a = y/\/6, 
h\ = h = z/^f?, and /i2 — 1 the eliminations 
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-3 (a (6 72 6 + 4 76 - 4 73 6^ - 4 6^ + 31 73 6 _ 64) 

_ 24 ^5 + 4 ^^ - 674 z^ + 1437 z"^ + 1355 z - 5292 
^ ~ 2(20^3 + 47^2- 117;z-604) 

generate the tenth-degree polynomial in z possessing the four real roots. These 
results are summarized in Table 5. 

3 Perturbations 

Any family of phenomenological potentials may be supposed approximated by an 
asymptotically harmonic Pade approximant (||) of a suitable degree t. Up to the 
K—i\i order a consistency of its subsequent perturbative treatment is guaranteed 
whenever the remainder is kept sufficiently small, v^ (phenomenological) j^^-j _ ^(Pado) j^^-j _ 

0{\^^) = 0{\^^J^). The stability of approximation requires that the more or less 
random poles in V^^^^^\x) are under &m control. In the case of our class of forces 
(H) this is most easily achieved by their unique |]16l partial-fraction re- arrangement 



v{(3,x) = cjv + y] y" — ^ + y] y] • (15) 

m=l j=l i'^ + emX'^y n=l k=l [i + ifn - '2gn)x'^ + g^X^f 

An instructive illustration is offered by the t = 2 example (|^). The maximal ad- 
missible range of perturbation of its couplings must be restricted by the condition 
of positivity of the denominator. This means that we must have / = /(A) > or 
/(A) = and (7 = ^(A) < or > /(A) > 4^(A) for A G (A„,in, A^ax)- 

We shall expand the perturbed wave functions in the same modified oscillator 
basis as above, 

oo oo t 

^{\,r) = B{r) J2 KiX) {r\n) = J2 ^n(A) (r|S„), i3(r) = J2 B^x^"- (16) 

71=0 n=0 d=0 



This leads to the infinite-dimensional Schrodinger equation 

[n{\) - E{X)V{X)]h{X) = 0. (17) 

It degenerates back to the finite-dimensional problems of preceding section in the 
unperturbed limit /3 — > (3^^^ . 

3.1 Phenomenological appeal of Pade oscillators 

The degree t in eq. (0) specifies also a half-bandwidth of our (2t + 1)— diagonal quasi- 
Hamiltonians [H(A) - E{X)V{X)]. This t = t[J(l), . . . , J(Mi), K{1), ..., K{N2)] 
grows rather quickly with all its arguments. Vice versa, the very first t's already 
offer a rich variety of possible shapes of the phenomenological potential (cf. Figures 
2 and 3). 

In the simplest t = 1 example (0) the numerically calculated A— or /3— dependence 
of energies exhibits a roughly quadratic shape, /? ~ {E + constY + const. For the four 
lowest states this is illustrated by Figures 1 and 4. The exceptional exact energies 
(marked by crosses) are scattered all over the coupling-energy plane. These quasi- 
harmonic points may be inter-connected by auxiliary lines (cf. Figures 5 and 6 and 
ref. ^7j). In a way, these lines generalize the harmonic spectrum to P^^^ ^ 0. Their 



distinguished feature seems to be an asymptotically almost equidistant and almost 
linear shape, resembling strongly their harmonic predecessor. 

Only the first few energies exhibit in fact a pronounced non-equidistant spacing. 
The onset of the almost equidistant behaviour moves only slowly up with the growth 
of P^^\ The approximate linearity of dependence of the n— th energy level En on the 
value of the coupling (3 is remarkable. We may expect that the first-order pertur- 
bation formulae will reproduce the t = 1 energies E{X) with decent precision in a 
broad interval of their A— dependence. 

It is well known that in the context of studies of double wells one of the big 
challenges to perturbation theory is posed by the related approximate degeneracy 
between the even and odd states. An explicit illustration of this phenomenon is 
provided by Table 1. It indicates that the long- lasting puzzle of perturbations of the 



quasi-degenerate spectra in the deep double wells may find one of its very natural 
resolutions in the present language since our formalism treats the states with differ- 
ent parity as perturbations of different systems. For example, in between the first 
two couplings of the Table [i.e., for (3 G (64.89. . . ,81.88. . .)] one should calculate 
the ground state energies as perturbations of E^^'' = 17 while the very close first 
excitations should be perturbations of E^^^ = 19 in another potential. In low orders 
the split of energies will probably remain disguised by errors of their separate per- 
turbative determinations but the related quasi-degenerate eigenfunctions themselves 
become clearly distinguished by their parity. 

3.2 Rayleigh-Schrodinger expansions near (3^^^ 

The measure A = /? — P^^'' of deviation of our perturbed Schrodinger equation from 



its zero-order form should be sufficiently small in the Kato's sense |]T8|. Then, any 
analytic A— dependence of the matrices 

v{x) = v{o) + Ai?(i) + A2r'(2) + . . . ^ ^(^) = ^(0) + A7i:(i) + a27i:(2) + _ _ _ ^ (ig) 

may be expected to imply the validity of the energy series @ and of its wave function 
counterpart 

h^ = hj{X) = hj{0) + Xhf + X^hf + ... . (19) 

This transforms our A— dependent Schrodinger equation into a set of its separate 
C(A'') components. At A; = 0, the unperturbed problem [7^(0) - E{Q)V{Q)]h{Q) = 
of preceding section is re-obtained. Next we get its 0{X) descendant 

[n{o) - E(°)r'(o)] h^^^ = [E(°)p(^) - n^^^] h{o) + E^^^v{o)h{o) . (20) 

In compact notation with abbreviations r(°) = [E^^^-D^^) -n'^^^]h{0), pf°> = V{0)h{0) 
and A4 = 7i(0) — E{0)V{0) this equation shares its form with all the subsequent 
0{X^) equations 

>l/^W=f('=-i)+EWpf°> (21) 

requiring only the further abbreviation 

fKD = [^(o)2)(2) + ^(i)p(i) _ H(2)]/I(0) + [E(o)d(i) + E(^)D(0) - n^'^]h'-'^ 
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and, in general, 
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As long as detA^ = each particular solution h^''^ may contain an arbitrarily large 
admixture of the zero-order column vector h{0). This is the well known renormaliza- 
tion freedom of perturbative wave functions in quantum mechanics. We get rid of it 
by the normalization /i^'^) = in each perturbation order A; > 0. Such a convention 
differs from the standard textbook recommendations but serves the same purpose 
and makes the solutions of our key eq. (0) well defined. 

3.3 Example: tridiagonal A4 

Let us choose t = 1, fix the nodal count q and accept an exact solution at /5 = (3^^^ as 
our illustrative ^°-' —superscripted zero-order approximation. The perturbed couplings 
P = P^^^ + A and B = B^'^^ + \ B^^^ + . . . with a small measure of perturbation A 7^ 
enter the infinite dimensional Schrodinger equation (^) or ([T7|). After appropriate 
insertions we get the t = 1 set of equations (pT]), 
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Only the first q + 2 components of fr-^^ = P^^-* h^^^ are nonzero since 
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df"^ = /^°^/3n, 
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In contrast, the compressed previous-order (i.e., already known) corrections r) 
only terminate in the first order, at /c = 1. 
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3.3.1 Ground-state illustration 



In the k = 1 and g = exemplification of equation (|2c 
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(25) 



the first two energy- dependent rows decouple from the rest. For the even-parity 
a = —1 they read 

As long as jS^^"* = 6 and Z?'-^-* = 1 their solution reproduces the current textbook 
first-order overlap formula for the energy, 



^(1) ^ /^exp(-x^)(l + x^)rfx ^ _6_ 
/^^exp(-x2)(l + x2)2rfx 11 ■ 

This test demonstrates the user-friendliness of our non-Hermitian recipe. 
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3.3.2 The first excitation 



In our preceding illustration we did not mark the cut-off M in A^. For the next, 
g = 1 state with M = 6 (chosen small for paedagogical purposes) we have to solve 
the set 
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The upper part of this equation (separated by the inset hues) stays decoupled. We 
can underhne that in contrast to the current textbook recipe we do not need any 
left eigenvector of Ai and still can define the k—th energy correction via the finite, 
{q + t + 1)— dimensional matrix inversion 
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In the second, lower part of eq. (p7D the cut-off M should grow to infinity in prin- 
ciple. The matrix-inversion evaluation of the remaining wave function components 
hg_lf_^_l, hg_l^_^_2, ... is much more difficult and acquires a purely numerical character. 
In a schematic representation h^'^^ ~ TZt^''~^^ + . . . the "unperturbed propagator" 7?- 
is a general, fully non-diagonal matrix. 



3.3.3 The higher excitations 

The size of the upper part of eq. (|23| ) grows with q (cf. ( |25| ) and (0)). Only its 
last, (g + 1)— subscripted line remains trivial. As long as Cg^i = 0, dg^i = and 
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g^^' > 0, it degenerates to an explicit definition of the {q + 1)— st coefficient 



in terms of the not yet specified energy. 
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Thus, in practice, our equation (p3| ) decays in the two equally difficult subsystems 
at the larger q. Its upper q + 1 rows are, fortunately, not as complicated as they 
look. Firstly, an immediate elimination of energies may help for g ^ 1. It proves 
unexpectedly easy since the array p^^^ is a left eigenvector of our non-Hermitian 
quasi-Hamiltonian 7W. The left action of this vector on eq. (^TJ) eliminates all h^^^^s 
and defines the k—th energy correction at any t, 
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With this knowledge, equation (|29| ) determines /ig+i. At the exact energy (pO]) an 
abbreviated right-had-side vector f^~^^ = rj^^^^ + E^'^^ p^\ m = 0, 1, . . . enters the 
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equation ( ]23|) whose first g + 1 separate rows become linearly dependent. We omit the 
very first one as redundant. Simultaneously, our normalization h'^^^ = annihilates 
a column in A4 and we get the reduced equation 



^4°) af) 



d'; 



(0) 



C2 a2 



V 



4°^ 



r(0) 






J 



'0 

,(fc) 



'''q-2 

V Vi / 



/ ^(^-1) \ 



(k-l) 

'2 



T^q-1 

. f{k~l) 
\ 1 



This gives a non-diagonal, upper triangular generalization 






'-0 

Ak) 



uik) 
'''q~2 



h 



(k) 



( 



1/c 



(0) 



7^l2 
1/4°) 7^ 



23 



7^2, 



V 






/ dk-^) \ 
'1 



T^q-1 
fik^-^) 

q 



(31) 



(32) 



of the current rule based on a diagonal unperturbed propagator. Unfortunately, any 
extension of the latter trick fails. For all t > 1 a different approach is needed. 



4 New recipe using triangular propagators 

Previous examples clarified an exceptional role of the (g + t)— th row in M. where 



-9+t 



d, 



q+t 



0. As a key source of simplifications at t = 1 this is not fully 



transferable to t > 1 [§ . At the higher t's, one is simply expected to pre-diagonalize 
the matrix M. in textbook spirit. This is a purely numerical step of course. In what 
follows we are going to describe a more analytic approach. It will lie somewhere 
in between the fully analytic (schematically, t = 0) and purely numerical (i.e., pre- 
diagonalization) extremes: Our unperturbed propagators TZ will be constructed as 
sparse, triangular matrices. 

The presentation of this material will be split in three parts. Firstly, using just 
t = 1 for simplicity, subsection W^ explains the idea of constructing TZ in the upper 
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triangular form. Secondly, subsection ^4.2| employs the next, t = 2 example and ex- 
plains the alternative approach using the lower triangular unperturbed propagators. 
We believe that this gives a clear guide to the general t's. Still, an abstract 
and detailed description of our innovated perturbation theory (with any t) is offered 
in the Appendix. The reason is that, building on the reader's experience with our 
previous examples, we may introduce a less transparent but much more compact 
notation. Moreover, we also relax there the immediate connection of our technique 
to some peculiarities (e.g., one-dimensional nature) of our illustrative example (^. 



4.1 Upper triangular propagators: t = 1 example 



We have split our eq. (^) into the separate rule 
solution (^2]) and the lower part 

7(0) \ 



), the upper part ( pTD with 



/ JO) 
"13+2 

^(0) 
Cg+3 



d 



a 



q+2 

(0) ,(0) 

g+3 "g+3 



.(0) 



,(0) 



^(j+4 "(7+4 



d 



(0) 
Af-1 



v 



JO) 
-M 



(0) 

«a/ / 



"■9+2 



h 



{k) 
q+3 



V.(k) 
'''q+2 



uik) 



( Jk-l)^ \ 



T, 



q+2 
~{k~l) 



T. 



<?+3 

~{k~l) 



T, 



5+4 



. Jk-1) 



(33) 



/ 



truncated at certain M ^ q + 2 and containing, on its right-hand side, a re-defined. 



~{k-l)* _ ~{k-l) 



-superscripted known quantity rg+2 = 'Tq+2 



JO) Jk) 
^q+2 "'q+1- 



For a clear explanation of our main idea let us now drop the superscripts and 
choose Q' = 1 and M = 6 again. This returns us back to our paedagogical example 
(p^) re-written now in a re-partitioned, equivalent square-matrix form 



/ 



C3 



V 



03 (is 
C4 (34 U4 
C5 05 (is 

C7 



\ 


(hA 




( ~ \ 

T-3 




h. 




7^4 




h. 


= 


^5 




h 




T^e 


/ 


[KJ 




^ ciK J 



(34) 



where /i2 is already known and a trivial last row c-jh^ = c-jh^ has been added. The 
trick is that we may now remove the cut-off completely. The infinite-dimensional 
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left-hand-side matrix 



/ .(0) 



Z 



-1 



-13+2 ^q+2 



(0) 



"g+2 



JO) (0) 

^(0) 



"g+3 

"9+4 



V 



(35) 



/ 



is regular. It may be inverted in an algebraic, non-numerical and cut-off-independent 
manner, "forgetting" our use of the vectors with M = 6, i.e., f/ = c-jh^ and /le+j = 
and fy+j = for all j = 1, 2, . . .. We may conclude that once we know our upper 
triangular propagator matrix Z we may pre-multiply by it equation (|34D from the 
left and get the final wave function defined by the formula 



,(fc) 



,(fc-i) 



(fc) 



,(0) 






m 



l,2,...,M-g-l 



(36) 



at any cut-off M ^ g + 1. Both its components have just an elementary form 



f^q+1 




^<2+2 




^q+l 




^i,M-q 


(fc-1) 

/^g+2 


= z 


^g+3 


1 


^q+2 




Zi,M-q 


(fc-1) 
f^M-2 




(fe-1) 
^M-1 




^AI-2 




ZM~q-2,M-q 


(fc-1) 

y /^A/-i / 












I ZM-q-l,M-q y 



(37) 



This is a key point of our considerations. Starting from the first omitted index 
m = M — q our equation (0) is an identity. At the first admitted (and exceptional) 
index m = 1 this equation defines, paradoxically, the right-hand side quantity Hm 
itself. Indeed, the pertaining left-hand side value hq+i is already known. 

In eq. ( p7|) an ascending recurrent evaluations may be recommended as giving, 

(fc-1) _ ^(fc-l) / (0) 



step by step, /i]v^_i' 



r 



M 



/c\l etc. Our recipe is complete. We may summarize: In 



a deeply anharmonic double well regime, the partial solvability of our unperturbed 
system and the reducibility of its propagator to a triangular matrix implies the 
feasibility of an innovated perturbation construction with propagator TZ of an upper 
triangular matrix form. 
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An additional, marginal remark is due. If needed, the recurrences (^) may be 
solved in a closed form 



.ly 



(fc-i) 

/^M-m-l (0) (0) (0) 

^M^M-1 ■ ■ ■ ^M-m 



where 



/ .(0) 



F = det 



''M-m "'M-m 



(0) 







.(0) 



,(0) 



(0) 



^M-m+1 ^M-m+1 "'M-m+1 









V 

due to the Kramer's rule. 



c 



(0) ^(0) 



(k-l) \ 



r; 



(fc-1) 
' M-m+1 






T. 



M 



(38) 



(39) 



4.2 Lower triangular propagators: t = 2 example 



Let us consider the pentadiagonal version of eq. (pl|) 



Oo 


"o 


df 








... 


c? 


ai 


d? 


d? 





... 


cf 


4" 


02 


rfW 





... 





4^' 


C3 


as 


4" 


<r ... 








C4 


C4 


04 


4" ... 



V 



\ 


( 4'^ ] 




( (k-i) \ 
^0 




Po 




(H?) 




rf 




pr 






= 


r?-) 


+ i^W 


pi'' 

Ps 

Pi 


.f-) 


^4 


/ 


\ ■ ) 




\ ■ / 




\ '■ ) 



J2] 



.{0} 



JO} 



with, say, 5' = 1, i.e., C3 = (^3 =0 and p\ = p\ = . . . = Q. After the insertion of 
energy (|30|) we sum the two right hand side vectors in one, rj + E^'^^pj = y , 
and omit the 3— subscripted row (between lines). Simultaneously our normalization 
JiW = (in parenthesis) annihilates the second column in A^. 

Moving further the first t — 1 (i.e., one) plus one (exceptional, (g + t)— th) columns 
of A^ to the right hand side and dropping the redundant superscripts ^''^ and ^''~^^ 
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we get the equivalent equation 



Z 



-1 





ho 





\ ''■ / \ 
with the appropriate auxihary f^^_,_i 
(while fX/+3 = ^M+4 = • • • = 0), with pi 
triangular and infinite matrix 





(40) 



Ca/+i hM-1 + c^j^i Hm and r^j+2 



[2] L 



{1} 



^ {2} 

ao etc, p\ 



d\ etc and with the lower 



f*' 










4" 


<i? 








a2 


4" 


dPi 






C4 


ci" 


04 


4=' 







4^' 


C5 


4' 


S? 


\ • 











J 

This pentadiagonal matrix is easily invertible since its main diagonal is all non-zero. 
Indeed, by construction, n ^^ t + q'm. S^ = Bt [Sn — ^t+q\ (n|r^*|n + 1) . The action of 
Z upon eq. (^) from the left gives our final wave functions 

/ h. \ 



t-2 

h3 



ha 



z 



( h\ 






hn Z 



p!" 
pi" 



/l5 Z 



\ pS;' / 



{2} 
P2 

{2} 



\ pi? ) 



(41) 



The seemingly redundant last two rows form in fact a core of the whole construction: 
With the well known left-hand side values of /lAf+i = ^a/+2 = they must be read 
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as the two necessary linear algebraic equations needed to determine the two "input 
parameters" ho and h^. 

4.2.1 An illustration 

A straightforward transition to the general t's does not require any new ideas but 
merely an appropriate "shorthand" notation. Its explicit description may be found 
in the Appendix. A nontrivial advantage of its more abstract language is that one 
can immediately work with a much more general class of Hamiltonians, say, with 
their non-Hermiticity and asymmetry related to the phenomenological absorption, 
etc. At the same time, even our simple Pade oscillators may provide a number of 
useful applications. 

Flexibility of our fairly weak assumptions may lead to a few non-standard con- 
structions. Imagine just an interpolation between two zero-oder models. Thus, our 
previous t > 2 solvable model (^ may be re-interpreted as the new potential 

V^,,,0^.^,M±(i^^)i^^^±^±(l^^)il!!!l (42) 

also solvable aX 6 = 0. With a new "small" parameter 6 = 1 — X and the new 
couplings 



/i 



(^) = 14.941997 ... - 16 ^ 1.05800, z/^^^ = 4.959146 ...-2^ 2.95915 



this leads to a methodically appealing linear interpolation between the two solvable 
cases. 

This may serve as a test of our method. In Table 6 the first-order precision 
compares well with the purely numerical exact energies. In a way paralleling our 
above t = 1 test (|2^) the usual evaluation of the t = 2 overlap integrals is much more 
tedious of course. Although the integrators in MAPLE [|T^ still manage and offer 
their evaluation, the immediate algebraic solution of our linear algebraic three-by- 
three equations proves, definitely, much preferable. 
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5 Summary 

Our present recipe treats a phenomenological Schrodinger equation in three steps. In 
the prehminary one we choose the potential in its Pade (or perturbed Pade) asymp- 
totically harmonic representation (|^). In the next preparatory stage the rational 
potential with 2t + 1 free parameters is assigned a suitable solvable zero-order ap- 
proximant with as many as t + 1 free parameters. We pick up the parity (or angular 
momentum) I and degree q of the unperturbed wave function and determine, alge- 
braically, all the parameters which are constrained by the solvability. In the third 
step we finally apply our modified or innovated Rayleigh-Schrodinger perturbation 
algorithm. In each order we 

• compress our knowledge of preceding corrections (say, in an auxiliary vector f 
given by formula (0)) and, if needed, 

• define immediately the new energy (by eq. (^)); 

• construct another auxiliary array (in general, vector 9 defined by recurrences 
( ^81) in the Appendix); 



choose a cut-off M ^ 1; 

satisfy the M— dependent model-space constraints (i.e., t -|- 1 linear equations 
(|0|) + dg) in general); 



• evaluate all the missing components of the wave function corrections (their 
general form is given by eq. (^Bf) below) . 

As a comfortable methodical alternative to the current prescriptions our innovated 
procedure admits a non-diagonality of unperturbed Hamiltonians and avoids the 
necessity of their pre-diagonahzation at A = 0. Its computational efficiency stems 
from its consequently recurrent character. 

For practical puposes, it is promising that our recipe is reducible, basically, just 
to a single recurrence relation per each perturbation order. This lowers the common 
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danger of a possible undetected numerical loss of precision, further suppressed here 
(and especially in the context of non-linear algebra in the zero-order constructions) 
by the high-precision computer arithmetics and programming in MAPLE ||19| . 

Originally, our choice of the illustrative rational potentials (^ has been motivated, 
mathematically, by a comparative smallness of their short-range perturbations. A 
posteriori, numerical tests clarified their phenomenological appeal. The flexibility of 
their shapes proved paralleled by the "fairly dense" occurrence of their partially solv- 
able bound states. Indeed, their observables (e.g., energies) exhibit often an almost 
linear or quadratic coupling-dependence. In a way extending the t = 1 observations 
by Gallas fl^ this supports a very good precision of perturbative predictions over a 



major part of the coupling space. We may expect a facilitated tractability of poten- 
tials with less common (e.g., multiple well) shapes. The frequent physical need of 
analysis of large variations of realistic forces seems to have found here an adequate 
computational tool. 
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Table 1. Low-lying spectra in the four deepest solvable double wells (p. 



/3(o) 


q. 


energies 


81.88 


3 


18.999999996 19. 22.765764732 22.765764788 
26.526337990 26.526338492 30.281295324 30.281297900 


64.89 


3 


17. 17.000000131 20.733677525 20.733679219 
24.459570379 24.459582299 28.176801248 28.176862466 


52.05 


3 


15.301693677 15.301695784 18.999974785 19. 
22.686594466 22.686760593 26.359595371 26.360391029 


49.91 


2 


14.999996593 15. 18.690932465 18.690972685 
22.369232872 22.369494489 26.032568100 26.033806209 


39.12 


3 


13.356890687 13.356934267 17. 17.000474393 
20.621974574 20.624839279 24.215073151 24.227701473 
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Table 2. Complete list of the t = 1 roots P^^\ 



q 


parity 


/?W 


excitation 


^(0) 





even 


6. 


ground state 


5. 




odd 


10. 


ground state 


7. 


1 


even 


8.8768943744 


first 


9. 






17.123105626 


ground state 


9. 




odd 


12. 


first 


11. 






26. 


ground state 


11. 


2 


even 


11.490856174 


second 


13. 






19.556337712 


first 


13. 






36.952806114 


ground state 


13. 




odd 


13.874580313 


second 


15. 






28.206711029 


first 


15. 






49.918708658 


ground state 


15. 


3 


even 


13.816182739 


third 


17. 






22.170398699 


second 


17. 






39.118906994 


first 


17. 






64.894511568 


ground state 


17. 




odd 


15.630566921 


third 


19. 






30.443898070 


second 


19. 






52.049183356 


first 


19. 






81.876351653 


ground state 


19. 
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Table 3. s— wave roots fi^^^ and u^^^ for potentials (|^ with f = g = 1 and g < 3. 





couplings 




coefficients 




Q 


/i(0) 


z/(o) 


u(0) 
"■0 


h^r^ 


/^f 


hP 





16 


2 


1 











1 


14.9420 


4.95915 


-3.48195 


1 








2 


14.0340 


7.91914 


8.18810 


-3.79751 


1 





2 


64.3015 


-2.40610 


1.02613 


1.82256 


1 





3 


13.2932 


10.8801 


-15.7092 


9.19559 


-3.91211 


1 


3 


62.7170 


0.883427 


-1.93699 


-2.48533 - 


-0.0989786 


1 
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Table 4. M — >■ oo convergence of s— wave energies for couplings (|T^) in (|^. 



M 


low lying spectrum 





11.422198 


— 


— 


— 


— 


— 


— 


1 


10.672913 


15. 


— 


— 


— 


— 


— 


2 


10.945092 


15. 


16.2817 


— 


— 


— 


— 


3 


10.944852 


15. 


— 


— 


— 


— 


— 


4 


10.944169 


15. 


18.5146 


— 


— 


— 


— 


5 


10.943697 


15. 


18.2954 


20.8470 


— 


— 


— 


6 


10.943435 


15. 


18.1857 


20.2068 


24.0949 


— 


— 


7 


10.943317 


15. 


18.1222 


20.0166 


23.6752 


28.065 


— 


8 


10.943284 


15. 


18.0876 


19.9294 


23.5621 


27.5340 


34.6709 


9 


10.943292 


15. 


18.0717 


19.8889 


23.5096 


27.4129 


31.3782 


10 


10.943316 


15. 


18.0670 


19.8730 


23.4838 


27.3538 


31.2387 


11 


10.943343 


15. 


18.0682 


19.8704 


23.4725 


27.3221 


31.1673 


12 


10.943366 


15. 


18.0718 


19.8743 


23.4696 


27.3063 


31.1255 


13 


10.943383 


15. 


18.0761 


19.8806 


23.4711 


27.3002 


31.1020 


14 


10.943395 


15. 


18.0798 


19.8871 


23.4745 


27.2999 


31.0907 


15 


10.943402 


15. 


18.0828 


19.8928 


23.4785 


27.3027 


31.0874 


16 


10.943406 


15. 


18.0849 


19.8971 


23.4822 


27.3068 


31.0889 


17 


10.943408 


15. 


18.0863 


19.9002 


23.4851 


27.3110 


31.0929 


18 


10.943408 


15. 


18.0872 


19.9022 


23.4874 


27.3147 


31.0977 


19 


10.943408 


15. 


18.0876 


19.9034 


23.4889 


27.3177 


31.1023 


20 


10.943408 


15. 


18.0878 


19.9040 


23.4899 


27.3198 


31.1083 
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Table 5. t = 3 roots u^^\ v^^^ and w^^^ for even parity, q < 2 and potentials (|T 



auxiliary root 




couplings 




coefficients 


z 


^m 


^(0) 


^(0) 


u(0) 


M°) 


— 


30. 


0. 


12. 


1. 


0. 


2.0534020780474 


58.107 


-14.804 


19.797 


1.452 


1. 


-11.95601933076 


30.088 


0.5698 


15.691 


-8.454 


1. 


3.9105205273467 


93.642 


-33.229 


33.945 


1.091 


2.258 


2.1316238057633 


86.526 


29.646 


15.572 


0.578 


1.231 


-4.8208818863125 


58.716 


-15.324 


23.422 


-6.313 


-2.783 


-11.9379718900079 


30.248 


1.295 


19.381 


33.908 


-6.892 



Table 6. Double test: Interpolation (^) between the two solvable potentials. 



"large" coupling 


A = l 


S=l-X=l 




ground state 


ffist excitation 


unperturbed energy 

the ffist correction 

k = 1 approximation 


11.0000 
-0.0449 
10.9551 


15.0000 
-0.0985 
14.9015 


Runge Kutta value 


10.9434 


14.6332 
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Appendix. Perturbation construction without left 
eigenvectors 

The main mathematical starting point of all our previous considerations was the 
linear algebraic zero-order Schrodinger equation Ai{E^^^)h^^^ = with elements such 
that 

Mm,m+t+j = 0, m = 0, 1, . . . , j = 1, 2, . . . 

Mn,n+t = 4*' > 1/^ > 0, n = 0, 1, . . . , n^no = t + q, (43) 

Let us now forget about any one- dimensional interpretation or Pade-oscillator origin 
of these matrices which are, by assumption, non-Hermitian. In such a case, it is 
natural to suppose that the left eigenvector of Ai is infinite-dimensional and ceases 
to be available. This appendix will show that (and how) our formalism remains 
applicable even without this auxiliary array. 

A. Notation 



The mathematical importance of the upper-diagonal dominance (|43| ) in A4 lies in its 



relevance for recurrences. Up to the exceptional n = tiq, each (i.e., the n—th) row of 



our eq. (pT]) may be read as an explicit recurrent definition of its "leftmost" unknown 




hl^\ The exception is unpleasant and its naive *— superscripted regularization 

f^i I ^ ^ w ioT n = nn = t + q 

4*^* = <| . (44) 

otherwise 

(k 

would produce a wrong value of /i2i+g without the limiting transition D ^ oo. 
A slightly more sophisticated recipe must be used. For its general formulation, 
the energy E^''' and elements h^''^ will be split in a pair of arrays C and ^. The 
former vector will be finite, containing just the t + 1 "difficult" components: energy 
E'C^) = Co, initial values hj = —Q, j = 1, 2, . . . , t — 1 (let us now prefer /ig = for 
normalization) and the exempted /i2t+g = ~Ct- AH the remaining unknowns will lie 
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in the infinite-diniensional vector ^ such that C,j = h- 



(k) 



0,1, 



i ^ q (notice 



that j = q would double-count the exceptional /i2t+g)- 

For the time being, let us leave the undetermined component of the new vector 
^ free. The presence of this new temporary variable z = ^t+q opens the possibility 
of a straightforward elimination of the vectors ^ as functions of the finite number of 
unknowns in (. Indeed, with an index j out of the interval j = 1, 2, . . . , t — 1 we 
may abbreviate all the j— th (i.e., leftmost) columns of our zero-order matrix Ai as 
respective vectors p^^^ distinguished by the braced superscript. We move them all to 
the right-hand side of eq. (^). This will reduce the left-hand side matrix Ai (acting 
on a vector) to the mere triangular submatrix denoted as Z from now on (and acting 
on a subvector). After the above simple-minded regularization (H3) with 



"o 



z -^Z* 



t/r ' < 





[t]* 








\ 



£ 



[t-2] Jt~l] At]* 



d 



d' 







V 



/ 



we may re-write our fundamental eq. (^Tj) in the presence of variable z 



Tn 



Z* 



6 



V ■ / 



Ti 



(fc-1) 



Tr 



(fc-1) 



V 



+ E0 



/ 



V '■■ J 



(45) 



The last column p^*J' must coincide with the irregular, 2t-|-g— th column of 7V1 without 
asterisk. The new form of our equation carries more information than before. The 
two auxiliary parameters (viz., D and z) were introduced. It easily follows from our 
construction that at any regularizing D ^ 0, our new equation remains precisely 
equivalent to its predecessor (^) if and only if the redundant variable z vanishes, 
z = 0. We may summarize: With the necessary limiting transition z — > deferred 
to the very end of our considerations, the most important computational benefit of 
the regularization D ^ may be seen in a completely recurrent solvability of our 
infinite-dimensional eq. (P^). Its first row defines ^o = h\ etc. 
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B. Reduction to a model space 



As long as, by construction, det Z* ^ 0, equation (|45|) could be solved by an imme- 
diate matrix inversion at an arbitrary t, D and 2, ^ = {Z*Y {t^^^^'> + I]j=o f^^^]- 
This would specify the infinite-dimensional left-hand side vector ^ as a finite sum. 



/ \ 
Co 


= 




\ ■ ) 




I ■ / 






Vo 

V • J 



(46) 



Accepting this idea, we may even try to determine each of the right-hand side com- 
ponents separately. This would be a tremendous simplification of the algorithm since 
the last t + 1 individual recurrences are order- independent, 



^ ^r 







dt'^ df* 








V 



\ 


^ vl," ] 




( Pf \ 




4" 


= 


p!" 


/ 


\ '■■ 1 




I ■ / 



0,1, 



(47) 



They just "reparametrize" a part of our zero-order Hamiltonian. Only the very first 
definition will vary with the order k, 



I 



d 



W 







4*-'^ 4*^* 








V 
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Still, it depends neither on the unknown energy E^'^'^ nor on the unknown coefficients 
/{(fc)_ We may say that it "compactifies" the previous, known results E^''~^\ }v-^~^\ 
E'^^~'^\ . . ., and "compresses" them into a new input vector 9^''~^^ = (Z*)^ t^''^^\ 
It remains for us to find the values of the t + 1 unknown parameters (j . For this 
purpose we return to the "forgotten" truncation conditions 



^MA 



^M+2 



ihi+t — 0. 



(49) 
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In the present notation these requirements read 

oik-i) ^t ^{j} A _ n 
n{k-i) I ^t Ai) A _ n 

^M+2 "T 2^i=0 '/Af+2Si — U 

/)(fc-l) ^i {j} /■ _ n. 

^M+t "T 2^i=0 'lAI+t^j — U 

and suppress the variabihty of parameters (. By construction, all these equations 
still depend on our auxiliary and, generically, nonvanishing variable z. Vice versa, 
the validity condition z = may only be re-introduced as an additional, explicit 
requirement ^t+g = 0, i.e., as the (t + 1)— st equation 

otg'^+i:viiiQ=o. (51) 

i=o 
The concatenated system (^TJ) + (^) of t + 1 conditions is our final model-space-like 
formula. It defines all the t -|- 1 parameters collected in the array C. 

At the correct value of energy (^) our auxiliary regularization variable z becomes 
identically equal to zero and the redundant equation ( pT]) may be omitted. Our 
model-space-like boundary conditions then degenerate to the mere t equations (pO]). 

C. Illustrations 

In the real fixed-point arithmetics an accumulation of errors may make the numerical 
value of z (i.e., the right-hand side of eq. (^Tf)) still slightly different from zero. 
Nevertheless, as long as the components of the latter equation itself are all of the 
order 0{D^^) for large D ^ 1 (cf. eq. (|5^)) the errors accumulate in the product 



Dz rather than in the quantity z itself. The choice of a sufficiently large D —^ oo 
settles the problem of errors. 

C.l. t = 1 

After we return to t = 1 for illustration, the last line of eq. (^) will represent our 
z = constraint (0). Similarly, the last line of eq. ( P7D plays the role of the second 
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constraint (0). Once we define the auxiliary two-column matrix rj, 

/ n \ 



r7^°> = (Z*) 






{0} 



r*\ — 1 



ag+2(0) 

V c5+3(o) y 

the inverse matrix in the pair of equations (|5lD and (0) remains triangular, 
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{0} {1} 

Vm+1 Vm+1 



n(k-l) 
^q+1 

flik-l) 



M>g>0, t=l. (52) 



This means that in any order k and at an arbitrary termination q and/or cut-off 
M < oo the definition of energies E'^''^ remains finite-dimensional. 



C.2. t = 2 

Above, we have already chosen g = 1 and M = Q for illustration at t = 2. Here, this 
gives the triangular regularized submatrix Z* which has the ordinary square-matrix 
form 
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t = 2. 



J 

Its elements marked by the superscript* are in fact arbitrary but, wherever possible, 

we shall keep them strictly equal to their nonzero unmarked values. Such a con- 
vention (with (ig = (ig , (Iq = (ig and rfg = rfg ) makes our construction less 



cut-off-dependent. Only the single element ^3 
nate. 



[2], 



D ^ remains really indetermi- 
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C.3. D^O 



Whenever the value of D vanishes in a "correct" hmit D ^ 0, matrix Z drops its 
asterisk and equation ( pi]) acquires its non- recurrent, pseudo- inversion form, 
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(53) 



Step by step, it defines the "upper" coefficients h2 \ h^ and hi as functions of 
the energy and of another undetermined parameter hi . Similarly, the first "lower" 
nonzero coefficient /ig is specified as a quantity which depends on all the "upper" 
coefficients plus on a new parameter hi . At the end, the correct values of our three 
unknown variables should be fixed by the "forgotten" fourth row and by the two 



truncation conditions h 



(fc)* 



h 



(kr 



0. 



C.4. L)>0 

In the example (p^D, the regularization Z{D = 0) —>■ Z*{D ^ 0) leads to the 
equivalent eq. (|45| ) if and only if the contributions proportional to z vanish, z —>■ 0. 
The situation remains virtually unchanged when we admit a growth of q. We may 
invert, in partitioned notation. 
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(54) 
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(with C = —Z3 + D ^wu ) and define the three order-independent and infinite- 
dimensional vectors 



^^0}^ r/^i}, ^^} 
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(55) 



It is an easy hnear algebra to show that the D— dependence of these components of 
the wavefunction corrections (cf. eqs. (|45D and (^41) ) implies their linear z— dependence. 
We may choose D ^ 1 and get the small numerical span of 2; = 0{D~^), i.e., only 
a small spuriosity in our tentative 2; 7^ wavef unctions, h^ {z) = hj (0) + 0{z). In 
the limit D ^ 00 we just return to the formulae studied above. Nevertheless, also 
any time before such a limiting transition, equations ( |5T[ ) and ( [50| ) define our last 
unknown parameters ( in terms of the vectors (|5^), via the three-dimensional matrix 
inversion 
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M > g > 0. (56) 



The occurrence of a zero matrix element is a peculiarity of the scheme (cf. eq. (p3|)). 
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Figure captions 

Figure 1. Coupling P vs. energy E for the first four bound states in potential (|l]) 

with B = 1. 

Figure 2. Three potentials (D supporting the exact ground state at zero energy. 

Figure 3. Three potentials (^ supporting the exact first excited state at zero energy. 

Figure 4. An empirical fit Y{E) = a + bE + cE"^ of the couplings (3 = (3{E). The 

thin and thick crosses denote the respective numerical and non-numerical levels at 

various I3^^\ 

Figure 5. The seven lowest energies for the seven lowest barriers Z?*^''-* in (|l]). The 

auxiliary "Gallas" lines connect the symmetric (upper curve) and asymmetric (lower 

curve) exact levels with growing q. 

Figure 6. Several "Gallas" lines in a bigger part of the E — (3 plane. 
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Figure 1 




Figure 4 




Figure 5 
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